Charge and spin density modulations in semiconductor quantum wires 
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We investigate static charge and spin density modulation patterns along a ferromag- 
net/semiconductor single junction quantum wire in the presence of spin-orbit coupling. Coherent 
scattering theory is used to calculate the charge and spin densities in the ballistic regime. The 
observed oscillatory behavior is explained in terms of the symmetry of the charge and spin distribu- 
tions of eigenstates in the semiconductor quantum wire. Also, we discuss the condition that these 
charge and spin density oscillations can be observed experimentally. 
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I. INTRODUCTION 

Relativistic quantum mechanics has predicted a rapid 
oscillatory behavior of free electrons, called Zitterbe- 
wegung, which is due to the interference between the 
positive- and negative-energy components in the wave 
packet^ This peculiar oscillation happens because the 
velocity is not a constant of motion even in the absence 
of any potential and has a fluctuating component, while 
the momentum commutes with the Hamiltonian. Re- 
cently, from the analogy between the band structure of 
semiconductor quantum wells and relativistic electrons 
in vacuum, it has been pointed out that electrons in 
semiconductors undergo the same oscillatory motion^^ 
The essential mechanism causing such a motion is the 
spin-orbit (SO) coupling in the semiconductor. In a two- 
dimensional electron gas confined in a heterostructure 
quantum well, two SO coupling effects are usually taken 
into account: the Rashba 4 and Dresselhaus^ SO coupling 
effects described by the following expressions 

HR = —?-{cT x py-<J y p x ) and H D = — -(<r y Py-<rxPx), 
m m 

(1) 

respectively, where er = (a x ,a y , a z ) are the Pauli matri- 
ces. The strength of each SO coupling is measured in 
terms of characteristic wave vectors kn and kjj, respec- 
tively. The Rashba term Hr arises when the confining 
potential of the quantum well lacks inversion symmetry, 
while the Dresselhaus term Ti.r> is due to bulk inversion 
asymmetry. In the presence of SO coupling, the velocity 
of the electron is not a constant of motion: 
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The spin precession a z (t) due to the SO coupling eventu- 
ally leads to a fluctuating velocity and consequently an 
oscillating motion of the electron. 

While the oscillatory behavior of the electron is dy- 
namic in the original proposal of the Zitterbewegung, it 
is highly plausible that one can observe static patterns 
of charge density oscillation in some properly-structured 
semiconductor samples. 



In our paper we investigate static charge and spin den- 
sity oscillations along a ferromagnet/semiconductor sin- 
gle junction quantum wire in the ballistic limit. The 
charge and spin densities across the system are calculated 
using coherent scattering theory. The observed oscilla- 
tory behaviors are explained in terms of the symmetric 
or antisymmetric structure of charge and spin distribu- 
tions of eigenstates in the semiconductor quantum wire. 
Also, we discuss the condition that these charge and spin 
density oscillations can be observed experimentally. 

Our paper is organized as follows: In Sec. [H] we investi- 
gate the properties of eigenstates of semiconductor quan- 
tum wires, paying attention to their symmetric structures 
in charge and spin density distributions. The effect of 
SO coupling on the structures in the presence of a con- 
finement potential is examined in detail. Section II I II is 
devoted to the numerical calculation of charge and spin 
density modulations in the ferromagnet/semiconductor 
single junction system in the presence of SO coupling. 
The main results are summarized in Sec. IIVI 



II. PROPERTY OF EIGENSTATES IN 
SEMICONDUCTOR QUANTUM WIRES 



Previous works&i have shown that perturbation theory 
cannot correctly explain the effect of moderate or large 
SO coupling on energy levels and properties of eigenstates 
in semiconductor quantum wires. Instead, truncating 
the Hilbert space to the lowest bandsSil or tight-binding 
models^ were successfully used to investigate the role of 
SO coupling on the transport through quantum wires. 
Here, we have implemented an exact numerical method 
of solving the system with arbitrary strength of SO cou- 
pling. Before introducing the method and discussing the 
properties of the states obtained numerically, we briefly 
review the structure of the system Hamiltonian and its 
symmetry in terms of transverse modes. 
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A. Hamiltonian and Symmetry Properties 

We consider a quasi-one-dimensional system of elec- 
trons in the presence of SO coupling: 



(3) 



where r and p are two-dimensional position and momen- 
tum vectors and m is the effective mass of the electrons 
in the semiconductor. The electrons are confined in the 
y direction by an infinite square-well potential of width 
L, 



V(r) = 



(|y| < L/2) 
oo (\y\>L/2), 



(4) 



We assume that the SO Hamiltonian H S o consists of Hr, 
and Hd , see Eq. QJ . In some semiconductor heterostruc- 
tures (e.g., InAs quantum wells) H.r dominates^, and in 
others (e.g., GaAs quantum wells) TLd is comparable to 
Hr—. A series of experiments^ has demonstrated that 
the strength of the Rashba SO coupling can be tuned 
by external gate voltages. Note that our choice of the 
square-well confinement eliminates the possibility of SO 
coupling due to effective electric fields coming from the 
nonuniformness of the confining potential. This kind of 
SO coupling should be considered in case of parabolic 
confining potential^ 

It is instructive to rewrite the Hamiltonian in Eq. @ 
in second-quantized form to gain insight into the effect 
of SO coupling on the transverse modes. To this end, we 
define Ck x nu/c\. to be the annihilation/creation oper- 
ators of the nth transverse mode \k x ,n, (m) with a wave 
vector k x and a spin branch index fi = ± in the absence 
of SO couplings. It is convenient to choose the spin po- 
larization axis n = (cos (p, sin ip) to be along the effective 
magnetic field due to the SO coupling for waves propa- 
gating in the x-direction such that 
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(5) 



with if = arg[fc/) + ikji]H In terms of these operators, 
the Hamiltonian is expressed as 



, c fc x n/j 



n// 



nn' u[i' 



(6) 



This expression shows that SO coupling leads to two ef- 
fects. On one hand, it lifts the spin degeneracy by shifting 
the energies of the transverse modes such that 



^nn{kx) 
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(k x - fik so ) 2 + (-^-) 



(7) 



with fc so = y/k B + k 2 D and A so = S 2 i s 2 /2m. On the 
other hand, it mixes the transverse modes according to 



the amplitudes k ri 



given by 
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(n =/= n' mod 2) 
otherwise 
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As shown in Eq. (JHJ), the SO coupling mixes transverse 
modes with opposite parities and possibly opposite spins. 
Therefore, the eigenstates of the Hamiltonian in Eq. © 
are generally not symmetric or antisymmetric. 

Still, the charge and spin density pattern can be sym- 
metric or antisymmetric: For simplicity, we focus on the 
case when only one of two SO coupling terms exists, that 
is, k R ko — 0. In this case, for any symmetric confine- 
ment potential satisfying V(y) — V(—y), the Hamilto- 
nian commutes with the so-called "spin parity" operator 
P s = P y (h ■ cr), where P y is the inversion operator for 
the y-component such that P y ip(x,y) = ip(x, —y)i& The 
eigenstates of the system should thus also be eigenstates 
of the spin parity operator, having the eigenvalues ±1, 
which is evident from the fact that P 2 = 1. Denoting 
v f , fc x ns(r) to be the wave function of an eigenstate for 
wave vector k x , subband index n, and quantum num- 
ber of spin parity s = ±1, one can show that the spin 
density components <S^ ns (r) = *^ ns (r)(m • cr)^/ k ^ ns (r) 
for m 1 n should be antisymmetric with respect to 
y = 0,2^ while the charge density and the spin com- 
ponent S™ ns for rh = n are symmetric. For exam- 



Pie, S'(y) 



k , ns 



(-y)(n ■ a)a z {h ■ cr)* 



k x ns \ 



-y) 



ns (—y)' This (anti)symmetric properties for rh = z 
persist even in case k R ko ^ 0, when there is no simple 
operator like P s commuting with the Hamiltonian. How- 
ever, as \kfi\ and \ko\ are tuned to come close to each 
other, the antisymmetric components decrease in mag- 
nitude and even vanish when \kn\ = \kn\ since the off- 
diagonal terms in the matrix w that couple modes with 
opposite spins vanish. Interestingly, for \kn\ = |&d|, the 
spin operator n • cr itself commutes with the Hamiltonian 
and the SO coupling couples only the transverse modes 
with the same spin (see Eq. I|8b|l) so that a common spin 
quantization axis can be defined for all the eigenstates. 15 

The symmetry properties discussed up to now are not 
restricted to the square-well potential case. One can show 
that the same reasoning works in the case of a harmonic 
potential.— All that one should do is to redefine the am- 
plitude k nn i according to the confinement potential. For 
the symmetric potential satisfying V{— y) = V(y), the 
amplitude k nn i has nonzero values only when n and n! 
have opposite parities and reproduces the same symme- 
try properties as discussed above. Hence, the choice of 
confinement potential does not lead to any qualitatively 
different effect. Therefore, we focus on the square-well 
potential since it is convenient for numerical calculations. 
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FIG. 1: (Color online) (a) Spectra of a quantum wire confined 
by a square- well potential in the presence of Rashba spin-orbit 
(SO) coupling: UrL — 1 (upper) and 3 (lower). The dotted 
lines correspond to the modes in the absence of SO coupling. 
The energy is scaled by Ei = 2mL' A ' ^ ne ground-state energy 
at k x = in the absence of SO coupling, (b) Spin expectation 
value {a y } for four low-lying states at /crL = 3. Same colors 
as in (a) are used to identify the levels, (c) Charge (black), 
spin-y (red), and spin-z (blue) density profiles along the y- 
direction for the lowest subband (n = 1) with s = +1 (upper) 
and -1 (lower) at k R L = 3 and E/E l » 1.74 2 . 



B. Numerical Calculation of Eigenstates 

It is possible to obtain the eigenstates numerically by 
diagonalizing the Hamiltonian Eq. JSJ in a truncated 
basis of transverse modes. Since we need eigenstates 
at a given energy, not a given wave vector, we have 
adopted another numerical method to obtain the eigen- 
states: First, prepare eigenstates at a given energy in the 
absence of the confinement potential^ and then numer- 
ically find a set of wave vectors k x such that the linear 
superposition of four eigenstates sharing the same wave 
vector k x satisfies free boundary conditions at the infinite 
walls, that is, vanishes at y = ±L/2. This method is su- 
perior to the diagonalization of the Hamiltonian 7ig nd in 
two aspects; one can obtain numerically exact eigenstates 
in this way, and an infinite set of evanescent waves, that 
is, eigenstates with complex k x at the targeted energy 
can be found systematically. Evanescent waves turned 
out to be important in transport through (semi-)finitc 
samples^ 

First, we consider the case where only Rashba SO cou- 
pling is present (hp — 0). Figure ^ shows the spectra of 
this system and the properties of the eigenstates for weak 
and strong SO coupling. Small SO coupling (kjjL < 1) 
splits the spin degeneracy according to Eq. 0) (see the 
upper figure in Fig. 0a)) and each spin branch s = ±1 
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FIG. 2: (Color online) Charge (black), spin- a; (green), spin-y 
(red), and spin-z (blue) density profiles along the y-direction 
for the lowest subband (n = 1) with s = +1 (upper) and — 1 
(lower) at (/crL, koL) = (a) (3,f ) and (b) (3, 2) and at energy 
EjE x w 1.74 2 . 



keeps its spin polarization direction approximately such 
that ((n • a-)) kxna = (vy) ksons ~ s as long as k x is not 
so large [k x < (2/crL 2 ) -1 ). Due to the antisymmetry 
of the spin distribution, the spin expectation values (<7 X ) 
and (a z ) strictly vanish irrespective of the SO coupling- 
strength. However, the spin density component ns (r) 
is finite and has opposite sign for different branches: 

s tns( r ) ~ - S k x n-s( r )- As k R L increases, on the other 
hand, the mixing of the transverse modes with oppo- 
site spins becomes more apparent such that level cross- 
ings happen even in the same subband for n > 1 (see 
the lower figure in Fig. 0a)) and the spin polarization is 
changed and even reversed compared to the case of weak 
SO coupling. Figure ^b) shows that at strong SO cou- 
pling (kuL — 3) and large k x (> (2knL 2 )~ 1 ) the spins 
of two branches in the lowest subband (n = 1) are po- 
larized parallel to each other: (<J y ) -, k ~ sgn(k x ) for 
both s = ±li This is due to mixing with higher-subband 
states having opposite spin, which manifests itself well in 
Fig. njc) where the charge density profile of the state 
s = — 1 of the lowest subband exposes the contribution 
from the higher subband n — 2. The charge and spin den- 
sity profiles in Fig.^c) also confirm the (anti)symmetric 
structure proven based on the symmetry of the system. 

Dresselhaus SO coupling Hd, related to Hr by a uni- 
tary transformation U = (a x + o y )/\/2, gives rise to an 
energy spectrum and charge and spin density profiles 
identical to those obtained in the Rashba SO coupling 
case except for the rotation of spin axes such that x — > y, 
y — > x, and z — > — z and the substitution of kn by ko- 

In the presence of both Rashba and Dresselhaus SO 
coupling terms, the spin polarization of the eigenstates 
is not along the x- or j/-axes any longer, since there are 
nonzero spin expectation values (a x ) and (cr y ), and its 
direction depends on the ratio kn/k]j. Nevertheless, as 
shown in Fig. [21 the spin density profiles are found to 
have (anti)symmetric structures similar to the previous 
cases: S^(y) = S^(-y) and S*(y) = -S*(-y). Note 
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that both spin- a; and y components are symmetric. The 
antisymmetric component S z (y) is observed to decrease 
in magnitude as kn approaches and eventually van- 
ishing when kn = kp. This behavior, consistent with the 
argument in terms of coupling matrix elements of w in 
Eq. (|8b(l , indicates that whereas Rashba and Dresselhaus 
SO coupling can induce an antisymmetric distribution of 
S z (y) separately, the presence of both of them will lead 
to a reduction in S z (y). 

In summary, SO interaction in the presence of a sym- 
metric confining potential can result in antisymmetric 
structures in the spin distributions of individual eigen- 
states, especially in the spin-z component. This phe- 
nomenon manifests itself strongly when either Rashba or 
Dresselhaus SO coupling is present and vanishes as both 
coupling strengths become equal to each other. 

III. CHARGE AND SPIN MODULATIONS IN 
THE PRESENCE OF SPIN-ORBIT COUPLING 

In this section we investigate the charge and spin den- 
sity modulations in a ferromagnet/semiconductor single 
junction quantum wire. Static charge and spin density 
oscillations are induced by injecting a spin-polarized cur- 
rent from ferromagnet to semiconductor. We make use 
of the coherent scattering theory to calculate the charge 
and spin density modulations in the ballistic regime and 
uncover the fact that the (anti)symmetric structure of 
eigenstates in semiconductor quantum wires studied in 
the previous section is deeply related to such density os- 
cillations. 



A. Scattering Theory 

The single junction system consists of the ferromagnet 
(F; x<0) and the semiconductor (S; a;>0) regions con- 
fined by the square-well potential of width L. While the 
electrons in the semiconductor region are governed by 
Eq. J2Jl, the Hamiltonian for the ferromagnet region is 

H F = ^ + V(r)-(h/2)(a z -l). (9) 

The spin-splitting energy h is assumed to be very large 
compared to other energy scales such that only the ma- 
jority spin state (spin up state) is available at the energy 
of interest. Here we have assumed that both regions have 
identical effective masses and lower band edges, because 
the impedance mismatch due to a difference in them does 
not affect the density oscillation in a qualitative way. In 
order to set up the scattering theory one should construct 
a set of eigenstates in both regions at a given energy 

F/S ' L-F / S F/S 

E. Denote (r) = e m * x $ m (y) to be eigenstates 
in the ferromagnet /semiconductor regions propagating in 
the x (e = +) or —a; (e — — ) direction at energy E, where 
m — (n, s) is a composite index with subband index n and 
spin index s and <& m (y) is a spinor for the state. Note 



that the wave vector k mx can be complex-valued in the 
case of evanescent waves. 

The boundary conditions at the interface (a: = 0) are 
then specified by the following two equations: 

^(r)U=o = * S (r)U=o (10a) 
v^ F (r)\ x=0 = v s x * s (r)\ x=0 (10b) 

for all y <G [— L/2,L/2}. Here the wave functions in 
both regions can be expanded in terms of the eigenstates 
as follows: tt» = E ro ^(r) i^cL^W 
and ty s (r) = J2 m c "^+m,( r ) w i tn complex coefficients 
c w">*) _ ^he velocity operators, v F = p x /m and v x = 
(p x — h{kji<7y + kD<7 x ))/m differ from each other due to 
the presence of SO coupling in the semiconductor. 

The SO coupling deforms the transverse modes such 
that they do not match in the ferromagnet and semi- 
conductor regions. It is thus impossible to reduce the 
scattering problem to a one-dimensional one. Instead 
we set up an infinite number of coupled linear equations 
for the coefficients c < - l ' r '*' by multiplying <&„(?/) on both 
sides of Eqs. IjlOaf) and (|10b(l and integrating them over y. 
The coupled linear equations are then solved numerically 
for the reflection/transmission coefficients at given 

(i) 

incident coefficients c m . In our study we focus on the 
injection of the electron through the lowest transverse 

(i) 

mode such that c; , N = 1 and the other incident 

(n=l,s=+) 

coefficients are zero. Since the contribution of evanes- 
cent waves in high subbands is very small in this case, 
the coupled equations can be truncated to be of finite di- 
mensions with negligible errors. By using the coefficients 
obtained in this way, we calculate and investigate the 
charge and spin density profiles, p(r) = | , I' F /' s (r)| 2 and 
5™(r) = * F/5t (r)(rh • a)^ F / s (r) in the next sections. 

B. Charge Oscillation due to Spin-orbit Coupling 

First, we investigate the charge density oscillation due 
to Rashba SO coupling and focus on the energy levels for 
which only two propagating states exist in the semicon- 
ductor region. Figure |3] shows the charge density profiles 
near the ferromagnet/semiconductor junction for various 
strengths of the Rashba SO coupling and clearly disclose 
the static oscillatory patterns perpendicular to the prop- 
agating direction. The patterns oscillate around the cen- 
ter line of the wire and the oscillation amplitude increases 
with the coupling strength k^lL. For large SO coupling 
(kjiL > 1.5) the electron does not stay any longer near 
the center line, forming high-density islands off-center. 
This oscillation is related to the fact that two propa- 
gating waves are eigenstates of the spin parity operator 
with opposite eigenvalues ±1. From this property one 
can prove that the charge density p(r) can be divided 
into symmetric and antisymmetric parts with respect to 
y = 0: 

p(x,y) = p s {y) + p a (x,y). (11) 
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FIG. 3: (Color online) Charge density profiles along 
the quantum wire in the neighborhood of the ferromag- 
net /semiconductor interface. A spin-polarized current is in- 
jected in the lowest transverse mode from the ferromagnet 
(x < 0) to the semiconductor (x > 0). The distances are 
scaled by L. Here we have set E/Ei ~ 1.42 2 and varied the 
strength of the Rashba SO coupling from kaL = 1 to 4. Ir- 
regular patterns near the interface are caused by evanescent 
waves whose contributions increase with k^. 
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FIG. 4: (Color online) (a) Period and (b) amplitude of the 
charge oscillation as functions of Rashba SO coupling strength 
knL for various energies E/Ei ~ 1.01 2 (red), 1.42 2 (green), 
and 1.74 2 (blue). The black solid lines are obtained from 
perturbation theory such that (a) the line goes like A = n/kR 
and (b) they show a parabolic dependence on ka. 



While the symmetric component p s (y) — p s (—y) is uni- 
form along the wire, the antisymmetric part p a (x,y) — 
—p a {x,—y) oscillates along the x-direction with period 
A = 2tt / Afc and generates the density modulations shown 
in Fig. [3J Here, Ak is the wave vector difference of two 
propagating modes. 

The oscillation period A depends on both the SO cou- 
pling strength and the energy of the incident electron, 
as shown in Fig.^Ja). At weak SO couplings the period 
is determined by the characteristic wave vector kn such 
that it decreases like A » ir/kR regardless of the energy. 
As kft increases further, on the other hand, the period 
increases with kji, having a minimum. In addition, the 
period gets longer at higher energies. The non-monotonic 
behavior of the period is another manifestation of mixing 
between transverse modes due to the SO coupling. Fig- 




x(L) 



FIG. 5: (Color online) Charge density profiles at kaL — 2 
and E/Ei « 6.68 2 . The density scale is the same as in Fig. [31 



ure^Jb) shows the dependence of the oscillation ampli- 
tude 2/ m ax on the SO coupling strength and the energy. 
Here the oscillation amplitude is given by the distance 
between the center line (y = 0) and the points where the 
density reaches its maximum. The parabolic dependence 
of i/max on kn at small SO couplings indicates that the 
charge density oscillation is a second-order effect of the 
SO coupling. As ka is increased, the amplitude does not 
increase any longer and instead saturates to y max ~ L/5 
as long as the condition on the number of propagating 
waves is not violated. 

At higher energies more than two propagating waves 
are involved in the scattering process in the semiconduc- 
tor region and the charge density oscillation pattern due 
to their interference cannot be specified by one period. 
Instead, some beating patterns appear as shown in Fig. [5] 
While the SO coupling still leads to oscillatory motion 
of the charge, oscillations of short and long periods are 
mixed and the pattern becomes complicated. 

Our observation is related to the prediction of Schlie- 
mann and Loss^ about the oscillatory behavior of free 
electrons in the presence of Rashba SO coupling: the cen- 
ter of an electron wave packet oscillates in the direction 
perpendicular to its group velocity v with the frequency 
2&r|v|, giving rise to the oscillation period ir/kn along 
the propagating direction. Controlled steady injection of 
wave packets with the same v can produce the static den- 
sity modulations we have observed. They also studied the 
dynamic oscillation of an electron with a given momen- 
tum k x in quantum wires, where the electron is prepared 
initially in the lowest-lying spin-polarized states. In this 
case, since the initial state is not an energy eigenstate, 
dynamic beating patterns due to the interference of waves 
with different energies are observed. On the other hand, 
our static density oscillation patterns are due to the in- 
terference of the energy eigenstates with a given energy. 
For that reason, such beating patterns and the resonance 
phenomena, leading to the sharp peak in the oscillation 
amplitude, are absent in our observation. 

We now include the effect of the Dresselhaus SO cou- 
pling on the charge density modulation. The unitary 
relation between two SO terms H.r and Tin discussed 
in Sec. [H] guarantees that the Dresselhaus SO coupling 
results in exactly the same charge density oscillation as 
the Rashba SO coupling. In the presence of both coupling- 
terms, however, the oscillation vanishes as the strengths 
of two terms become equal to each other. Figure 
shows that the amplitude of the oscillation decreases 



6 



0.08 - 

3 

a 0.04 - 
E 

j ^ , ^ ► 

05 1 

ko/kR 

FIG. 6: (Color online) Charge oscillation amplitude as a func- 
tion of the ratio kn/kn for kn,L = 1 and EjE\ ~ 1.74 2 . 

with kjj/kfi for finite kn and the oscillation disappears 
when kn = fc^. This tendency is quite similar to the 
dependence of the antisymmetric distribution of S z (y) 
on the SO coupling. It also supports the close relation 
between the charge density oscillation and the antisym- 
metric structure of eigenstates in semiconductor quan- 
tum wires. 

Before closing this section we address the experimen- 
tal observability of the charge density modulations. Since 
the oscillation patterns depend on the energy of the inci- 
dent electron, it should be checked whether the superpo- 
sition of charge density patterns from different energies 
can weaken the overall density oscillation. Figure 0Ja) 
shows that for large SO coupling {HrL > 1.5) the period 
strongly depends on the energy. Therefore, in the case 
of a large voltage drop across the junction the injection 
of electrons with a wide range of energy will smear the 
density modulation by superposing density patterns os- 
cillating with different periods. For small SO coupling, 
on the other hand, the phase as well as the oscillation pe- 
riod pattern are weakly dependent on the energy, leading 
to the superposition of commensurate density patterns. 
In addition, the states below the fermi level are irrele- 
vant to the charge density oscillation because the levels 
are filled incoherently such that their charge distributions 
are always symmetric with respect to y = 0. Hence the 
condition for the observation of charge density oscillation 
is that either the SO coupling strength or the applied 
voltage across the junction should be small enough. Un- 
der this condition, charge oscillations can be detected via 
high-resolution scanning probe microscopy techniques 16 
by imaging the charge density directly or by measuring 
the change of conductance along the quantum wire as a 
function of the position of a tip that can deplete charges 
under it^ 



C. Spin Oscillation due to Spin-orbit Coupling 

In the presence of SO coupling the electron spin pro- 
cesses in the effective magnetic field that depends on the 
momentum. Spin density patterns thus reflect such spin 
precession as well as the charge density distribution dis- 
cussed in the previous section. Figure shows spin den- 
sity profiles along the quantum wire if only Rashba SO 
coupling is present. The spin precessing patterns in 5 z (r) 
indicate that the effective magnetic field is approximately 
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FIG. 7: (Color online) Spin density distributions, S y (r) (up- 
per) and S z (r) (lower) for knL — (a) 1 and (b) 2 in the 
same configuration as used in Fig. |3] Here we have set 
EjEx « 1.74 2 . Note that S*(r) = 0. 
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FIG. 8: (Color online) Spin density distributions, 
S x (r) (upper), S y (r) (middle), and S"(r) (lower) in the same 
configuration as used in Fig. |31 Same parameters and color 
scaling as in Fig.0are used except that knL — 2 and knL — 1. 



perpendicular to the electron's moving direction. The 
spin distribution 5 z (r) is not symmetric with respect to 
y = in the semiconductor region and this behavior gets 
stronger as the SO coupling increases. It is because the 
spin pattern follows the charge density oscillation that 
is also asymmetric. While S x (r) = for all r, S y (r) 
is nonvanishing inside the semiconductor, increasing in 
magnitude with the SO coupling strength. Unlike S z (r), 
S y (r) get biased to positive values as ku is increased, in- 
dicating a spin polarization along the y-direction. This is 
a direct consequence of the spin polarization of the lowest 
modes at large k x and k^ (see Sec. Ill B|) . 

The Dresselhaus SO coupling also gives rise to similar 
spin profiles; the spin axes x and y are exchanged, and 
the spin distributions are mirror-reflected with respect 
to y = 0. In the presence of both coupling effects, on 
the other hand, the spin density patterns, with all the 
components nonvanishing, become distorted as shown in 
Fig- El Slanted patterns appear due to the absence of any 
symmetry with respect to y = 0. As long as &r ^ 
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the spin is preferentially polarized along either x- (when 
fco > kji) or y- (when k,R < ko) direction. For kn = kr> 
an unbiased spin precession takes place for all the spin 
components. 

Under the same condition found for the detection of 
charge density modulations, the spin density distribu- 
tions can be experimentally observed. Spin-polarized 
scanning tunneling microscopy^ with magnet-coated 
tips or optical techniques such as the spatially re- 
solved Faraday rotation spectroscopy^ can produce 
high-resolution images of spin density profiles. 

IV. CONCLUSION 

We have investigated the symmetry properties of eigen- 
states in semiconductor quantum wires and observed 
the charge and spin density modulations along quantum 
wires consisting of a ferromagnet/semiconductor junction 
by using the coherent scattering theory. We have shown 
that the Rashba or Dresselhaus SO coupling terms induce 
charge density oscillations perpendicular to the propa- 



gating direction and that the (anti)symmetric structure 
of charge and spin distributions in eigenstates is deeply 
related to the oscillations. Charge and spin density oscil- 
lations can be experimentally observed as long as the SO 
coupling strength or the voltage drop across the junction 
are small enough. 

Impurities, dirty interfaces, and electron-electron in- 
teractions that are not included in our study affect the 
transport in low-dimensional systems. Although these ef- 
fects will therefore influence the form of the charge den- 
sity modulations and make it more homogeneous, the 
charge density modulation may still be observed in suffi- 
ciently clean and low-density samples. 
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